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Abstract 

Different holographic dark-energy models are studied from a unifying point of view. We compare 
models for which the Hubble scale, the future event horizon or a quantity proportional to the Ricci 
scale are taken as the infrared cutoff length. We demonstrate that the mere definition of the 
holographic dark-energy density generally implies an interaction with the dark-matter component. 
We discuss the relation between the equation-of-state parameter and the energy density ratio of 
both components for each of the choices, as well as the possibility of non-interacting and scaling 
solutions. Parameter estimations for all three cutoff options are performed with the help of a 
Bayesian statistical analysis, using data from supernovae type la and the history of the Hubble 
parameter. The ACDM model is the clear winner of the analysis. According to the Bayesian 
Information Criterion [BIG), all holographic models should be considered as ruled out, since the 
difference ABIG to the corresponding ACDM value is > 10. According to the Akaike Information 
Criterion (AIG), however, we find AAIG < 2 for models with Hubble-scale and Ricci-scale cutoffs, 
indicating, that they may still be competitive. As we show for the example of the Ricci-scale case, 
also the use of certain priors, reducing the number of free parameters to that of the ACDM model, 
may result in a competitive holographic model. 
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I. INTRODUCTION 

Observations from Supernova (SN la) []]], large scale structure cosmic microwave 
background [sl, the integrated Sachs- Wolfe effect |4|, baryonic acoustic oscillations ^ and 
from gravitational lensing {g] suggest that our current universe is in a state of accelerated 
expansion. In most investigations the origin of the acceleration is attributed to a mysterious 
component with a negative pressure, called dark energy (DE). The preferred candidate for 
this entity is a cosmological constant A. The favored cosmological model is the A-cold- 
dark-matter (ACDM) model which also plays the role of a reference model for alternative 
approaches to the DE problem. Alternative DE models have been developed since the ACDM 
model, although it fits most observational data rather well, suffers from two main and lasting 
shortcomings, namely: the low value of the vacuum energy and the coincidence problem j^]. 
In order to address these two problems, the cosmological constant has been replaced by a time 
varying quantity, resulting in dynamical DE models. The best studied models here are scalar 
field models which comprehend, e.g., quintessence js], K-essence jsl and tachyon fields 10|. 
In most investigations DM and DE are considered as independent substances. However, there 
is no reason to neglect interactions in the dark sector. Including the possibility of a coupling 
between DM and DE is the more general case and gives rise to a richer dynamics [11]. So 
far, DE and dark matter (DM) manifest themselves only through their gravitational action. 
This circumstance is the motivation for unified models of the cosmological substratum in 
which one single component plays the role of DM and DE simultaneously. A Chaplygin 
gas represents the prototyije of a unified model 3], but also bulk- viscous models belong to 
this category (see, e.g., [13(| and references therein). Among the many different approaches 
to describe the dark cosmological sector, so called holographic DE models have received 
considerable attention 
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16| . According to the holographic principle, the number of degrees 



of freedom in a bounded system should be finite and related to the area of its boundary [Ij 
Based on this principle, a field theoretical relation between a short distance (ultraviolet) 
cutoff and a long distance (infrared) cutoff was established [15[. This relation ensures that 
the energy in a box of size L does not exceed the energy of a black hole of the same size. 
Applied to the dynamics of the Universe, L has to be a cosmological length scale. Different 
choices of this cutoff scale result in different DE models. If one identifies L with the Hubble 
radius H~^, the resulting DE density, corresponding to the ultra-violet cutoff, will be close 
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to the observed effective cosmological constant |15|. The possibihties of the_^particle and 
the event horizons as the IR cutoff length were subsequently discussed by Li [16], who found 
that apparently only a future event horizon cutoff can give a viable DE model. However, 
afterwards it was recognized, that a cutoff at the Hubble scale may well result in a realistic 
cosmological dynamics provided only, that an interaction in the dark sector is admitted 18 1. 
More recently, yet another possibility was proposed, namely a cutoff scale, given by the Ricci 



scalar curvature 



19 



20| . resulting in so-called holographic Ricci DE models. Holo grap hic 



DE model have been tested and constrained by various astronomical observations 2l|, in 
some cases also including spatial curvature {22 1. 

A sp ecial class are models in which holographic DE is allowed to interact with DM 
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23|-|26| . In this article, we compare three different approaches to (generally interacting) 



DE models. We study the relation between the energy density ratio of DM and DE and 
the equation-of-state (EoS) parameter for general interactions in each of these cases. Also 
solutions for the interaction free limits and scaling solutions (constant energy density ratio) 
are obtained as special cases. We clarify the role of potential interactions in the dark 
sector. In particular, we point out that any interaction model, including also the non- 
interacting limits, introduces relations between the EoS parameter and the matter content, 
which constrain the admissible dynamics. All the cases are confronted with observational 
data from supernovae type la and from the history of the Hubble parameter. A statistical 
analysis is performed on the basis of both the Akaike Information Criterion {AIC) and the 
Bayesian Information Criterion [BIC). It turns out, that a cosmological constant remains 
the favored DE candidate. 

This paper is organized as follows: In Section |TT] we provide the basic relations for an 
arbitrary cutoff scale. In Section IIIII we consider the Hubble-scale cutoff. The cutoff at the 
future event horizon is studied in Section llVt including an analysis of the non-interacting 
limit and the possibility of scaling solution. The same features are investigated for the 
Ricci-scale cutoff in Section |Vl In Section |Vl] we present the results of a Bayesian statistical 
analysis for each of the cases. Finally, a summary of the paper is given in Section I VI II 
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II. BASIC EQUATIONS 

We assume the cosmic substratum to be described by a pressureless matter component 
with energy density pm and a holographic dark energy component with energy density p^- 
The Friedmann equation for the spatially flat case then is 



?,H^ = %'nG{pra + PH) 



(1) 



Both components are admitted to interact according to 



Pm + ^Hpm = Q 



(2) 



and 



PH + 3Hil + w)pH = -Q. (3) 

Here we have introduced the equation-of-state (EoS) parameter w = where pn is the 
pressure associated with the holographic component. The total energy density p = Pm + Ph 
is conserved. The Hubble rate changes as 



H 



-HU1 + 



w 



1 + r 



dhiH 
d\n.a 



w[a) 



1 + r(a) 



where r = ^ is the ratio of the energy densities of both components and 



w 



1 + r 



w 



(4) 



(5) 



is the total effective EoS of the cosmic medium. The djTiamics of the ratio r is determined 

by 



Pm Ph 
^Pm Ph , 

Introducing here the balances ([2]) and ([3]), we arrive at 



?>Hr (1 + r) 



Following 
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w ^ Q 

1 + r 3Hpm_ 



16| we shall write the holographic energy density as 



PH 



L2 



(6) 



(7) 



(8) 



The quantity L is the infrared (IR) cutoff scale and Mp = l/VSvrG is the reduced Planck 
mass. The numerical constant determines the degree of saturation of the condition 



L^PH<MliL, 



(9) 



which is at the heart of any holographic dark energy model and which states that the 
energy in a box of size L should not exceed the energy of a black hole of the same size. To 
establish a specific model, the length scale L has to be specified. Applied to the dynamics 
of the Universe, L has to be a cosmological length scale. Three different choices of L have 
been considered in the literature: the Hubble scale, the future event horizon and a scale 
proportional to the inverse square root of the Ricci scalar. Each of these choices will be 
discussed separately in the subsequent sections. All relations of the present section are 
generally valid. 

Differentiating the expression ^ and inserting the result into the balance ^ yields 

T = 9_ = 2^-?,H{l + w) , (10) 

PH ^ 

where F denotes the (generally time-varying) rate by which ph chang result of the 

interaction. From the outset, there is no reason for Q to vanish. Putting Q = establishes 
a specific relationship between w and the ratio of the rates ^ and H. Any non-vanishing Q 
will modify this relationship. We shall investigate here the general dynamics and afterwards 
consider the non-interacting limits as special cases for each of the choices for L. 
With Q from (fTOj) the dynamics ([7]) of the energy density ratio r becomes 



-3H (1 + r) 



w 2 L 
1 + 



r 3HL 

The interaction free limit is characterized by (cf. eq. ([7])) 



(11) 



Q = =^ r = ^— -3h] =3Hrw (12) 
l + r \ L I 

with a generally time-dependent w. Different choices of the cutoff scale L result in different 
expressions for the total effective EoS parameter w^-^-^ in ([5]) and in different relations between 
w and r. It will turn out that for a Hubble-scale cutoff the ratio r is necessarily const ant, while 
any, generally time varying, w ^ requires Q ^ 0. For both the other choices, the future 
event horizon and the Ricci-scale cutoffs, there exist relations between w and r. In particular, 
in both these cases a constant ratio r requires a constant EoS parameter w. In the following 
sections the three choices for L are treated separately. 



5 



III. HUBBLE-SCALE CUTOFF 



General relations 



For L = H ^ the holographic DE density is 



PH = 3c^M^H\ (13) 



Differentiating this expression and applying Eq. (jl]) yields 



Consequently, 



PH = -3H{1 + ^ j . (14) 

Ph + 3H{1 + w)ph = S^^Y+V ■ ^^^^ 
This means, the source term Q, given by the right-hand side of the balance equation ([3]), is 
fixed: 

Q = -3Hw'ffpm ^ ^'^"" = -W^- (16) 
This constitutes a direct relation between the EoS parameter and the so far unspecified 
interaction quantity Q. An interaction with Q > is essential to have a negative EoS 
parameter. Without interaction there exists only the solution with w = 0, i.e., for the 
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interaction-free case this cannot be a dark-energy model 

Combining (JTj) and ( JT6l) (or ([5])), we find that necessarily r = const. Any interacting 
holographic dark-energy model based on (fT3l) is characterized by a constant ratio of the 
energy densities. We emphasize that a constant value of r does not imply that the EoS 
parameter is constant as well. The ratio remains unfixed. It can be freely chosen. 

Combining the Friedmann equation ([1]) and the expression f|T3|) . we obtain the following 
relation between r and the saturation parameter c^: 



^2 



1-C 



,2 



l + r) = l ^ r = — . (17) 



c 



A specific value of r fixes the saturation parameter. For r ^ ^, i.e., 25% DM and 75% DE, 
we find ~ |. For the (generally time varying) deceleration parameter one has 

■ ' " Ul-^). (18) 



H'^ 2 V Hpr, 

In the interaction-free case we recover the Einstein - de Sitter value q = \- The condition for 
accelerated expansion is Q > Hpm- To describe a transition from decelerated to accelerated 
expansion, Q has to change from Q < Hpm to Q > Hpm- 



B. Variable EoS parameter 



The freedom in the choice of can be used to estabhsh a viable cosmological model. 

I) that the interaction rate T is proportional to a power of the Hubble 

rate, equivalent to 

The quantity /i is an interaction constant. Different interaction rates are characterized by 
different values of n. A growth of the parameter is obtained for n > 0. In the spatially 
fiat background the ansatz (|T9|) corresponds to a total equation of state parameter 



^-'"-'^--^(i-V = -^(^Y"^' ■ (20) 



P \HoJ \Po, 

The energy density (for n ^ 0) is that of a generalized Chaplygin gas [27j |: 

p = po [/i + (1 - /i) a-'-/'] ^ H = Ho[f. + {l-f.) a-'-/'] . (21) 



Since from ffTSl) and (fT9l) it follows that 



p = i (1 - 2go) , (22) 



where go is the present value of the deceleration parameter, the Hubble rate can also be 
written as 

^ = (0 [1 - 2go + 2 (1 + go) a'^-l^] . (23) 

For the special case n = 2, this expression is similar to that for the ACDM model. The 
results of a statistical analysis are given in subsection IVIBt encoded in figures |5] and [61 



C. Constant EoS parameter 

The simplest dynamics is obtained for a constant EoS parameter, corresponding to n = 
in ( 12U]) . With Q = 3fj,Hpm we have = — p = const. Under this condition the balance 
([3]) takes the form 

PH + SH{l-fi)pH = 0. (24) 

To avoid an increase of pn, i-e., a phantom behavior, we have to require p < 1. This 
guarantees = — > —1 for the effective EoS. The deceleration parameter ( !T8|) 
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reduces to g = | (1 — 3/i). The condition for accelerated is g < =^ 3/i > 1. Exponential 

expansion is obtained for g = — 1 =^ fi = 1. From equation flMj) the solution for the 
densities is 

~ P„^ ~ P ~ a-3(^-^) (25) 

with 

2 2 , , 

a ~ and = — — . (26) 

3(l-p)t ^ ^ 

It is well known, that any power-law solution for the scale factor corresponds to a scalar- 
field representation for the cosmic substratum with an effective exponential potential V. In 
the present context we have 

V = Voexp[TX{<p-<l)o)] , (27) 

with 

A = ^/2A^^G~(lTry(^T^I^ . (28) 

For A to be real, w"^^^ > — 1 is required. This coincides with the condition for avoiding a 
phantom-type behavior, mentioned below (!24|) . 

IV. FUTURE EVENT HORIZON CUTOFF 
A. General relations 

Since accelerated expansion cannot be accounted for by holographic DE with an Hubble- 
scale cutoff without interaction (cf. (|T8l)). Li [16*] replaced the Hubble scale by the future 



event horizon Re- This approach has received considerable attention subsequently ^ 
With L = Re the holographic DE density (E]) is 

PH = —BT^ , (29) 

where 

^ . X ^ . r dt' r da' 

is the future event horizon. Differentiating (129|) yields 

pH = -2^ph ■ (31) 



R 



E 
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From (130!) we obtain 



Re Re 



Combining relations f l3ip and f l32p results in 

Ph = -2{H- R^') ph . 

For the left-hand side of the balance we obtain 

Ph + 3H{1 + w)ph = [(1 + 3w) H + 2R-'] pn • 

Comparing this equation with ([3]), we may read off 

2 



Q = -H 



1 + 3u; + 



ReH 



Ph w = ~- 



1 + 



+ 



H 



(32) 



(33) 



(34) 



(35) 



Here we have used that Re = c\/ (1 -\- r)/ H, a relation that allows us to eliminate c and to 
write 

/ 1 -u T r R2_f7'2 -| 

(36) 



1 + r 

ReH = ReoHqJ— ^ r = tq + (1 + ro) 



1 + ro 



RlH' 



1 



Notice that tq = y3q\' where ^2^0 = • To have an energy transfer from DE to DM 

one has to require 

(37) 



Q > 4^ m;<--(1 



ReH 



The dynamics of the ratio r is determined by ([7]). With (!35|) in ([7]) we may also write 



r = —Hr 



l + 2{REHr--^ 



or 



-H{l + r) 



w 2 
1 + 3- + 



1 + r ReH_ 

One expects the ratio r to decay during the cosmic evolution, i.e., 



1 / 2 

r<0 ^ u7>--(l + r) 1 + — — - 
3 \ ReH 



(38) 



(39) 



(40) 



This amounts to a lower bound for the EoS parameter. Combination of the upper bound 
fl37|) with f l40|) provides us with following range for the total effective EoS: 



w 

< < 



1 + 



(41) 



3V ReH J 1 + r 3 (1 + r) V ReH ^ 
Notice that q = 1(1 + 3y^). For a given w{a), Eq. ( 139|) determines the ratio r(a). A 
constant r implies that w is constant as well. This property is substantially different from 
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the case of the previous section where a cutoff at the Hubble scale with a constant ratio r 
still left room for an arbitrarily time- dependent EoS. The case of a constant r for the present 
future-event-horizon cutoff will be considered in subsection IIVDI below. 
The dark-energy balance (E]) an be written as 

PH + SHil + wffj)pH = (42) 

with an effective EoS 

where the superscript E denotes "event horizon". This quantity does not directly depend 
on w. However, the ratio r that enters ReH is determined by w via Eq. ( I39ll . Notice also, 
that this effective EoS for the DE component is different from the total effective EoS of the 
cosmic medium which is In the previous Hubble-scale-cutoff case both these quantities 
were identical. Via (143|) . a constant wfj-r necessarily implies a constant r and vice versa. 



B. Scaling energy-density ratio 



To solve the dynamics of the present model, we shall assume a power-law dependence for 
the energy-density ratio [29 1 

r = roa"^ . (44) 
Under such circumstances we have r = —C,Hr which, inserted into fl5^ . provides us with 



w 



1 + r - + -y/T+r 

c 



(45) 



The interaction rate that corresponds to a dynamics with 

T = Hr 



IS 



(46) 



C y/1 + r_ 

The first inequality in ( HTl) is satisfied for any ^ > 0, while the second inequality in ( HTll . 
equivalent to Q > 0, requires 



(47) 



ca/1 + r 

For any ^ > 1 this implies an upper limit for the ratio r. The ansatz (lHj) was introduced to 



approach the coincidence problem 



As was argued in 



29|, any value ^ < 3 makes the 
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coincidence problem less severe than in the ACDM model for which one has ^ = 3. In the 

following we consider the cases ^ = 1, 2, 3 separately. 

(i) Integration of with ^ = 1 leads to a holographic energy density 

4/c 



Ph = Pma 



(^ = 1) 



(48) 



With pr, 



rpH the total energy density becomes p = (1 + r)pH, equivalent to a Hubble rate 

2/c 



H = Hoa 



-2 Vro + a 



The free parameters are Hq, tq — 



(e = l) 



(49) 



l-f^o' 



and c. 



(ii) The corresponding relations for ^ = 2 are 



Ph = Pma 



Vro + 1 + 1 , 



2/c 



(e = 2) 



and 



H = Hoa 



yJrQ + a? + a 



. l/c 



(e = 2) 



(50) 



(51) 



(iii) For = 3 one has 



pH = Pma 



-2 



A/ro + + vc 
V'^o + 1 + 1 



1 4/(3c) 



(e = 3) 



and 



H = Hoa 



V^o + a^ 



i/ro + + 



2/{3c) 



(e = 3) 



(52) 



(53) 



With the observationally preferred values for and tq (see table [Tll in section |VT] below), the 
condition (HTll is satisfied for all the three cases, i.e., the energy transfer at the present time 
is always from DE to DM. But these cases are not only quantitatively but also qualitatively 
different. The case = 1 has the best-fit value = 1.14. Together with the preferred value 
Vq = 0.25 this corresponds to a present effective EoS parameter wfjrj in (H3|) of wfjjr{a = 
1) ^ —0.89. Even in the limit a ^ 1 the parameter wfj^j will remain larger than —1 and 
both the energy density pn in ( HHl) and the Hubble rate in (H9l) will decay in the long-time 
limit. The situation is different for ^ = 2 with the preferred values = 0.73 and tq = 0.51. 



Although the present value of the effective EoS parameter w 



effy 



-0.97 is larger 



than —1, it will cross the phantom line for a > 1. Both the energy density in eq. ( 150|1 and 
the Hubble rate in (|5T|) increase with a. For = 3 with = 0.42 and Tq = 0.95, the EoS 
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parameter (jl3l) has a present value of w^jf{a = 1) ^ —1.07, i.e., it is of the phantom type 
today and will remain so in the future. Also here, both the energy density pu in fl52|) and 
the Hubble rate fl53p are increasing functions of the scale factor. 

A statistical analysis provides us with figures [TJ El and [3 in subsection I VI CI below. Notice 
that all models with a cutoff at the future event horizon suffer from the serious drawback 
that they cannot describe a transition from decelerated to accelerated expansion. A future 
event horizon does not exist during the period of decelerated expansion. This drawback 
manifests itself here also in the inequality (H7j) which is clearly violated if r is larger than 
a certain threshold value. In other words, a matter dominated era with r ^ 1 cannot 
be recovered in this context. None of the Hubble rates (H^ . (|HT]) and ([S3D approaches a 
matter- type behavior H oc a"^/'^ for a <^ 1. 



C. Non-interacting limit 



Different from the previous Hubble-scale cutoff, there exists a non-interacting limit with 
accelerated expansion in the present case. In this special situation Eq. (135!) reduces to a 
relation w = w{r). With (5 = one has from (!35|l that 

2 „ „ 2 



ReH 



3w 



ReqHq = Cy/TT 



3wo 



(54) 



For Wo ~ — 1 and ro ~ | we have ~ | as in the Hubble-scale cutoff case. But notice that 
to derive this value for the Hubble-scale cutoff, neither an assumption Q = nor an EoS 
parameter wq ~ —1 was necessary. Together with ( 136|) we obtain r = r{w), 



1 + r 
1 + ro 



1 + Swo 
l + 3w 



r = ro + (1 + ro) 



1 + 3wo 
l + 3w 



1 



Alternatively, 



w 



1 + ro 



(55) 



(56) 



It is obvious, that for any Wq ~ —1, the parameter w remains always smaller than — |, 
demonstrating again the impossibility of a matter-dominated period in this context. This 
circumstance is visualized in a r — w plane in Fig. [T]for two values of Wq. The corresponding 
dependence r(a) is shown in Fig. [2] for the parametrization w = Wq + Wi{l — a), where Wq 



and 1^1 are constants 



3l|. 
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FIG. 1: Event-horizon cutoff, non-interacting limit. The plot shows the ratio r versus the parameter 
w for two values of wq with tq = 3/7. 

D. Scaling solution 



According to (138|) . a stationary solution r = tq requires 



1 + 



ReqHq 



(57) 



ifp^ 1 + ro 3 

Comparison with ( H3|l shows, that in this case the total effective EoS parameter coincides 
with the effective EoS parameter w^j- Notice, that r = is a separate requirement here. 
For the Hubble-scale cutoff a constant r was a property of the model. A further difference 
is the following. While for the Hubble-scale cutoff the constant ratio r was compatible with 
a time-dependent equation of state, the EoS parameter w is necessarily constant for the 
present configuration. It is completely determined by r and the saturation parameter c. A 
similar property holds for the ratio in fl57|) . This ratio is fixed as well. Recall, that 
for the Hubble-scale cutoff the corresponding quantity could be chosen freely. Notice that 
according to the last relation in ( 1571 ) the saturation parameter is given both by tq and wq. 
For the Hubble scale cutoff it was given by r = tq, independently of the EoS parameter. 
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1.0 



FIG. 2: Event-horizon cutoff, non-interacting limit. The plot shows the ratio r versus the scale 
factor a for two values of wq for ro = 3/7 and w\ = 0.31. 



For constant r the energy densities behave as 

pocpmOcpHOC a-2[i-{-RBoifo)-i] _ 

Integrating Priedmann's equation, we obtain the power-law solution 



a oc ii-i/(«i5o«o) , 



1 



1 



1 - 1/{ReoHo) t ■ 



The deceleration parameter becomes 



(58) 



(59) 



(60) 



ReoHo cVl + ro ' 

As expected, q is always negative. For holographic dark energy with an infrared cutoff at the 
future event horizon, a stationary solution for r is not compatible with a matter dominated 
phase. We recall that an event horizon does not exist for decelerated expansion, i.e., it is 
impossible to describe a transition from decelerated to accelerated expansion within a model 
with future event horizon cutoff. 
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The equivalent scalar field potential for this case coincides with ( 127|) and ( 128|1 where now 
w'^^-^ is determined by from fl57|) . 



V. RICCI SCALE CUTOFF 



General relations 



The role of a distance proportional to the Ricci scale as a causal connection scale for 



perturbations was noticed in 32|]. As a cutoff length in D 



33 



34 



36 



5 models it was first used in 19 |. 

The Ricci scalar is given 



Observational constraints were obtained, e.g., in 
by -R = 6 ^2iJ^ + if The corresponding cutoff-scale quantity is = -|, i.e. 



R 



PH = 3c^M}- = a[2H' + H) , (61) 
where a = Upon using (j4]) we obtain the expression 

for the holographic dark-energy density. It is remarkable, that the EoS parameter explicitly 
enters pn- With SH"^ = StcGph (1 + r) we find 

l = jil + r-3w) w = ^il + r)-^ . (63) 

Already at this stage it is obvious, that a constant value of w necessarily implies a constant 
r and vice versa. The time dependence of w is related to that of r by r = 3w. The last 
relation in flB^ may be used to express in terms of the present values of w and r: 

2 

— = l + ro-3wo, r = ro + 3{w -Wq) . (64) 

Recall that for the Hubble scale cutoff it was only the ratio r = tq which determined the 
saturation parameter. Here, c is again related both to tq and Wq, although by a relation 
that is different from the corresponding relation for the cutoff at the future event horizon. 
Differentiating (|6T|) yields 

Ph = a(4HH + . (65) 

From (jl]) one obtains 

2 V l + ^y 2 \l + r (l + r)V ^ ^ 
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Introducing (jl]) and (166|) in (l65l) . the latter becomes 



Ph = a 



-6H^ 1 + 



+ -if ^ 1 H 

1 + r/ 2 V i + r 



w 



wr 



1 + r {l + rY 



■ (67) 



From the definition fl6T]) with (j4]) we have 

3if (1 + w) = a 
Adding up the expressions (1671) and (l68l) we find 



6H'il + w)-lH'il+w) + ^ 
2 V 1 + r 



(68) 



Ph + 3H{1 + w) ph = -Q , 



where 



Q 



3H 



1 + r 



rw — 



w 
H 



pH ■ 



(69) 



(70) 



This general relation for the interaction term is a property of the model, determined by the 
ansatz fl6Tl) . No further assumption has entered here. The DE balance ([3]) may then be 
written as pn + 3H (l + w^j) Ph = with 



w 



eff 



w 
1 + 



w + 



H 



1 + To + 3 [W - Wq) 



(71) 



The superscript R stands for "Ricci" and tq 



i-n. 



— . Recall that the total effective EoS 



parameter is given by ([5]). Eq. (jlj) for the present case takes the form 

dlnH 3 1 + ro + 4 (w - fwo) 



(72) 



d\na 2 1 + ro + 3 {w — wq) 

For a given dependence w = w{a) the last equation provides us with H{a). The deceleration 
parameter q = —1 — ^ becomes 



Its present value is 



_ 1 1 + rp + 3 (2w - Wq) 
^ 21 + ro + 3(w- Wq) 

1 1 + ro + 3wo 



(73) 



(74) 



2 1 + ro 

The transition between phases of decelerated and accelerated expansion is given by g = 
in ( I73|) . Denoting the corresponding value of the EoS parameter by Wq we have 



1 



1 



g - - g (1 + ro) 



(75) 
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B. The CPL parametrization 



For a given dependence w{a), relation ( !75|) provides us with the value Og of the scale 



factor and, equivalently, the redshift at which the transition occurs. With the CPL 
parametrization w = Wq + Wi{l — a) we find 

Iwq 1 1 + ro 1 



3l| 



a. 



l + + ^. = --1. (76) 



' 2Wi ' 6 Wi ' ttq 



In order to have a transition before the present time, i.e. Og < 1, the condition 

2wi Q Wi 
has to be fulfilled. For t^i > this reduces to 

1 , 

Wo < -- (1 + ro) , 

which coincides with the condition for accelerated expansion at the present time (cf. 
Eq. fl7i|) ). For wi < we have the opposite inequality, so this case has to be excluded. 

1 the other SNIa datasets as well as the 



42|. 



With the exception of the Constitution dataset, a_ 
CMB and BAO data prefer positive values for wi 

Positive values of Wi are also expected from the behavior of the energy density ratio 
r = ro + 3{w — wq). For a < 1 one expects r > ro because the matter fraction is supposed 
to evolve towards smaller values with increasing a. This is guaranteed for w > Wq for a < 1. 
For negative w this means, that in the past w was less negative than the present value wq. 
With w — Wq = wi{l — a) this is realized for wi > 0. 

A further point is the direction of the energy transfer. From (1701) it is obvious that for any 
constant negative w the source term Q is positive, implying a transfer from DE to DM. For 
a time varying w the positivity of Q requires rw — aw' < 0. For the CLP parametrization 
this becomes rw + awi < 0. At the present time this reduces to Wi < —tqWo, providing us 
with an upper limit for Wi. 

With the CPL parametrization, eq. fl72|) can directly be integrated. The resulting Hubble- 
function is 



3 l + rQ + iUQ+Aini 

H = Hna ^ i+'-o+S"! 



1 + ro + 3wi (1 



_ 1 l + TQ-SufQ 

(77) 



1 + ro 

The free parameters are Hq, wq, wi and ro. No assumption about Q was made to obtain the 
expression (I77|) . The results of the statistical analysis are displayed in figure [TO] in the first 
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part of subsection IVIDI Notice that for the observationally preferred data the mentioned 
hmit wi < —tqWo seems to be violated. 



C. The interaction Q = 3H(3pH 

Relation (!70|) is valid for any interaction. In this subsection we combine it with the 
frequently used interaction model Q = SHPpH- Together with the second relation of f l64p 
this gives rise to a first-order differential equation for w which has the solution 

lu — s — (u + s) Aa^ 



w = 



6 



1- Aa' 



(78) 



where 



u 



ro - Swq + 3/3 



tq + 3wo + 3/3 



12/3 (1 + ro - 3wo) 



(79) 



and 



A = 



V 



V + s 



(80) 



With w{a) from (ITHI) . the ratio r{a) is explicitly known as well via In figures [3] and H] 
the functions w{a) and r(a), respectively, are shown for three values of /3. In both figures 
the (not included) non-interacting limits /3 = are almost indistinguishable from the cases 
P = 0.001. With w{a) and r(a) explicitly known, Eq. (jlj) for the Hubble rate may be 
integrated. The result is 



H 



a 2 



na — m 



n — m 



where 



m = 1 + ro - -{v - s) 



k = - {u- s) 





n 



1 + '^o - 2 + ■5) 



A 



I 



6 



{u + s)A . 



(81) 



(82) 



(83) 



The parameters here are ro, Wo, Hq and (5. The statistical analysis results in figure [TT] in 
the second part of subsection IVIDi Observationally preferred is a small, positive value of 
/3, i.e., a transfer from DE to DM. 
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FIG. 3: Ricci-scale cutoff. Tlie plot sfiows tlie parameter w versus tlie scale factor a for three 
values of the interaction parameter /3 for wq = —1 and ro = 3/7. 

D. Non- interacting limit 

While, according to (!70l) . a constant EoS parameter w has the non-interacting hmit w = 
for which behaves as dust, a non- interacting hmit (3 = does exist for w 0. This is 
different from the Hubble scale cutoff case, where the non-interacting limit has w = always. 
In the present case 

^ w dliaw 

Q = ^ rw = — r = — . (84 

H a ma 

In this limit the total effective EoS becomes 

^0 ^ _f.„_.,„„a ri , ^ ^ • (85) 



1 + r roa-^^'o-^wo) [i + _ 3^^^] _ 3^^^ • 

It is interesting to consider the limits of w{a) and r(a) for small and large values of the scale 
factor. At high redshifts we obtain 

w 

r -> ro - 3wo , -^-^ ^ ^ (a < 1) . (86) 
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FIG. 4: Ricci-scale cutoff. The plot shows the ratio r versus the scale factor a for three values of 
the interaction parameter (3 for wq = —1 and rg = 3/7. 

The far-future limits are 

w 1 

r ^ , — >wq- -ro (a > 1) . (87) 

1 + r 3 

The limit ( 186|) implies that the total cosmic medium behaves as dust, thus reproducing 
an early matter-dominated era. Assuming tentatively ro ~ | and Wq ~ —1, the ratio r 
approaches r y for a ^ 1. This value is only ten times larger than the present value tq. 
Recall that for the ACDM model there is a difference by about nine orders of magnitude 
between a value of r taken at the recombination epoch and its present value. In this sense, 
the coincidence problem is considerably alleviated for the the present model. On the other 
hand, in the far-future limit a ^ 1 the ratio r approaches zero as for the ACDM model. 
Apparently, the far-future EoS can be of the phantom type for wq — |ro < —1. The non- 
interacting limit /3 = of ( IHTj) becomes 

Ho Y 3wo - ro [1 + ro - 3wo\ 
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E. Scaling solution 



A scaling solution is realized for (cf. Eq. (^^] 

r = ^ yjR = 9_ , (89) 

"■^^ 1 + r 3Hp„ ^ ' 



With (1641) we then have li; = as well. Under this condition „^ = u = const. An 

^ ^ oil pm 

interaction with Q > is essential to have a negative EoS parameter. The second relation 
in flM|) coincides with ([5]) for the Hubble scale cutoff, only that there a constant ratio r does 
not necessarily imply that w is also constant. For the density ph we then have 

p^ = |i72(l + 3/i) . (90) 

The solution for the densities are identical with f l25|) and fl26|) of the Hubble scale cutoff. 
To avoid an increase of pu-, i-e., a phantom behavior, we have to require p < 1. This gives 
rise to the inequality (1 + ro) > |. On the other hand, the condition to have accelerated 
expansion is /i > |, equivalent to (1 + tq) < 1, i.e., the allowed range is 

^<c2(l + ro)<l. (91) 

If we combine a specific energy density ratio, say r = ro ~ | with an observationally 
suggested present EoS w ~ —1, we obtain ^ ^ ~ 0.46, i.e., an exact value which is inside 
the range (19T|) . In other words, with the observational input '"o ~ | and Wq ~ —1 the model 
is completely specified. In particular, the observational data for ro and wq determine the 
degree of saturation of the holographic condition ([9]). Notice that the value found here is 
smaller than the corresponding value for the Hubble-scale cutoff. Also for the present case 
the scalar field representation coincides with that for the Hubble-scale cutoff. In particular, 
we obtain the potential fl27|) with fl28|) . The dependence fl27|) with fl28l) is characteristic for 
all the stationary solutions. 

VI. STATISTICAL ANALYSIS AND OBSERVATIONAL CONSTRAINTS 
A. General description 

Each of the holographic models described in the previous sections provides us with a 
specific prediction for the behavior of the Hubble parameter as a function of the scale 
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factor. The relevant expressions are ( !23|) for the Hubble-scale cutoff, ( H9|) . ( 15T|) and (153|) 
for different cases of the event-horizon cutoff, flTTl) for the Ricci-scale cutoff using the CPL 
parametrization and (IHTj) for the Rici-scale cutoff with an interaction proportional to the 
DE density. We summarize them again as function of the redshift z = ^ — 1. 



Hubble radius. 



;i-2go) + 2(l + go)(l + ^: 



3n 
2 



(92) 



The free parameters are Ho, qq and n. In a. first step, the Hubble parameter Hq is de- 
termined by minimizing the three-dimensional function. The remaining parameters 
then are go and n, for which we perform a statistical analysis in subsection IVI B[ 



Future event horizon with ^ = 1. 

^3/2-l/c /l + ro(l + ^) 



H(z) = Ho(l + z) 



^/r^^Tz)TT+l_ 
Vro + 1 + 1 



2/c 



(93) 



Since ^ is fixed, we have Ho, ro and c as free parameters. Ho is obtained as in the 
previous case. The free-parameter space then consists of ro and c. The correspond- 
ing analysis is described in subsection IVI CI For the physical interpretation see the 
paragraph following Eq. f l53|) . 



Future event horizon with ^ = 2. 

,1-i/c ll + ro{l + z) 



H(z) = Ho(l + z) 



ro + 1 



^ro{l + zy + l + l 
Vro + 1 + 1 



l/c 



(94) 



As in the previous case, the free parameters are Hq, ro and c. See again subsection 
IVI CI and the paragraph following Eq. ( l53l) . 



Future event horizon ^ = 3. 



ro + 1 



Vro(l + ^)=^ + l + l 



2/(3c) 



(95) 



See subsection IVI CI and below Eq. f l53|) also here. 



Ricci scale with CPL parametrization. 



3 l+''o+"o+4t^l 
H{Z) = Ho{l + Z)^ l + ro+3-l 



l + ro + 3a;iY^' 
1 + ro 



1 l+'"0 — 3"0 
'2 l+ro+3t^i 



(96) 
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The free parameters of this model are Ho, ro, wq and wi. In this case, the minimum 
value of the four- dimensional x^-function is used to determine both Hq and tq. Then, 
a two-dimensional analysis is performed for wq and Wi, as described in subsubsection 
IVIDII 

Ricci scale with interaction Q = 3H(3pH. 



H{z) = Hoil + z) 



II 1- 



nil + z) 



3 Im — kn 
I 2 rans 



(97) 



n — m 

Here, one has Hq, tq, w and (3 as free parameters. We fix w = —1 and determine Hq 
along the lines already described for the previous cases. The statistical analysis for vq 
and (3 is the subject of subsubsection IVI D 21 

In order to probe the above models against observations, we consider four background 
;ests which are directly related to the behavior of the function H{z): the supernova type la 
it, the age of the very old galaxies leading to a direct measure of the H{z) function [43 1. 
the CMB shift parameter R 4J| and the baryonic acoustic oscillations BAO 45|. We shall 
present the results for a combined analysis of these four tests. 



The supernova-type-Ia test is based on the luminosity distance function 

c /"^ dz 



D, 



;i + ^)- 



(98) 



Ho Jo ^/H{z)' 

where c is the velocity of light. The observational relevant quantity is the moduli distance, 
given by 

/i = m-M = 51n(-^ )+25, (99) 



Mpc 

where m is the apparent magnitude and M is the absolute magnitude of a given supernova. 
In what follows we shall use the data set of the Union2 sample 47 1. 
The CMB shift parameter 



R 



dz 



Hiz) 



(100) 



depends on the matter density parameter VLmo and the redshift at the decoupling between 
matter and radiation, z^ = 1090. Observationally, R = 1.725 ±0.018 [48]. Baryonic acoustic 
oscillation (BAO) are characterized by the parameter 

x/n;;;;; ri n dz ^'/^ 



A 



[H{z,)]y^[zJ, H 



;ioi) 
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where Xb = 0.35. Observationally, A = 0.469 ± 0.017 Another test we will use is the age 
of the very old galaxies that have evolved passively. Our analysis is based on the 13 data 
for such objects listed in reference 49 |. 

For each of these observational tests we evaluate the fitting function x^, given by 

X' = Y1 2 ' ^ > (102) 

i=i ^« 

where ef^ stands for a theoretical estimation of the ith data of a given quantity (moduli 
distance, parameters R and A, H{z)), and stands for the corresponding observational 
data, (jj being the error bar. From this statistical parameters we can construct the probability 
distribution function (PDF), 

P{xf) = Ae-^'("')/2^ (103) 

where A is a normalization factor and x-^ denotes the set of free parameters of the model. 
The joint analysis is performed by summing up the different contributions, obtained for 
each of the tests separately. A one-dimensional PDF for a given parameter can be obtained 
by integrating over the remaining ones. All one-dimensional estimations will be made at 
2cr (95% confidence level), but we will also display the la, 2a and 3a contours for the 
two-dimensional PDFs. 

In order to gauge our results, we perform at first the analysis for the ACDM model. In 
the fiat case, there are two free parameters, Hq and Qmo (or, alternatively, ^ao). The best-fit 
value for Hq, based on the supernovae data, is Hq = 70.05 km/sec Mpc. For our purposes, 
the most important quantity is Qmo for which we find (at 2a) Qmo = 0.27lo'g2- In table [T] 
we display the values for Xmin G^ch model and for each observational test. The values 
for the parameters R and A are only shown for illustrative purposes. Here, the fitting is 
remarkable good due to the fact that there is just one observational point. The values in the 
last column of table [T] seem to indicate that all the models are competitive. But one should 
keep in mind that the ACDM model has just two free parameters, whereas the number of 
free parameters is three or four for the other models. In our analysis, however, we have fixed 
one or two of them, respectively. The last two columns of table IIIII show, how the statistical 
analysis penalizes the existence of additional parameters. In the following subsections we 
list the results for the different models in more detail. 
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SN 


H{z) 


CMB 


BAQ 


Total 


ACDM 


542.71 


8.07 


0.0 


1.02 X 10^20 


550.78 


Hubble(0„o = 0.25) 


542.59 


8.02 


3.27 X 10^15 


1.02 X 10-20 


550.61 


Hubble(QmO = 0.3) 


542.59 


8.02 


1.27 X 10-is 


1.61 X 10-15 


550.61 


Event horizon = 1) 


544.49 


8.53 


1.08 X 10"^^ 


5.08 X lO^i'^ 


553.03 


Event horizon = 2) 


542.79 


8.09 


5.29 X 10-^"^ 


4.82 X 10"^^ 


550.88 


Event horizon(^ = 3) 


542.95 


8.19 


5.26 X 10-1'^ 


3.09 X 10^20 


551.15 


Ricci (CPL param.) 


543.44 


8.01 


1.37 X 10"^^ 


2.02 X IQ-^^ 


551.45 


Ricci(Q = 3H(3ph) 


542.66 


8.02 


2.53 X 10"^^ 


1.70 X 10^20 


550.67 



TABLE I: Xmm values for each of the models and different tests. 



B. Hubble-scale cutoff 

In this case, there are three free parameters: Hq, qq and n (cf. Eq. (192!) ). The mass 
density has no direct connection with these parameters. This is important to stress since 
the CMB shift R and the BAO parameter A explicitly depend on the present mass density 
VLmo- The evaluation has been made for the cases f2mo = 0.25 and f2mo = 0.3. We minimize 
the three-dimensional function for a given set of observational data and use the value for 
Hq obtained in this way in order to reduce the space of parameters to two. This procedure 
is appropriate, since the PDF for Hq is sharply peaked around its maximum value. It will 
also be applied to the other cases to be studied in the following subsections. Generically, 
Hq = 70 km/sec Mpc for the supernovae data and Hq = 68 km/sec Mpc for the H{z) 
data. For Vlmo = 0.25 the results for the remaining parameters are go = ~0.61^q;o4 and 
n = 1.98;;^:^^, while they are go = -0-Q8-om and n = 2.47^[!;^^ for n^o = 0.3. The first 
case is very close to the ACDM model. The one- and two-dimensional PDFs for go and n 
are displayed in figured for Q^no = 0.25 and in figure [6] for Q^no = 0.3. Note the existence of 
remarkable oscillations in these PDFs. 
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FIG. 5: Hubble scale cutoff with rimo = 0.25. Left panel: one-dimensional PDF for qq. Center panel: 
one-dimensional PDF for n. Right panel: contour plots for the la, 2a and 3a confidence levels. 




FIG. 6: Hubble scale cutoff with ftmo — 0.3. Left panel: one-dimensional PDF for qq. Center panel: 
one-dimensional PDF for n. Right panel: Contour plots for the la, 2a and 3a confidence levels. 



C. Future event horizon cutoff 

The free parameters are the present values of the energy density ratio vq and the Hubble 
rate Hq, as well as the saturation parameter c and the power ^. We use the same procedure 
as in the previous case in order to determine Ho. Since we fix ^ to be either ^ = 1 or ,^ = 2 
or ^ = 3, we are left with two free parameters, tq and c (cf. Eqs. ( 1^ . (IMj) and (^5^). 
In table HIl the parameter estimations for tq and (the more convenient) are shown for 
different values of C,- There is a clear tendency: as ^ grows, the density ratio increases 
and the value for decreases. In figures [TJ [8] and [9] the one and two dimensional PDFs 
are exhibited. The physics of each of the models has been discussed in the part following 
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e = i 


e = 2 


e = 3 




"-'•^3-0.04 


'-'•^-'--o.or 


0.95l°i? 


C2 


1 14+0.05 
■'-•-'-^-0.05 


n 70+0.06 
^- '"^-0.06 


0.42^°:°^ 



TABLE II: Parameter estimations for vq and for different values of ^. 




FIG. 7: Future event horizon cutoff with ^ = 1. Left panel: one-dimensional PDF for tq. Center panel: 
one-dimensional PDF for c^. Right panel: Contour plots for the la, 2a and 3a confidence levels. 

Eq. (153|) . Although these models are unable to reproduce a matter-dominated phase at 
high redshifts, the results of the statistical analysis do not differ substantially from those of 
the models with a correct H oc (1-1-2;)^/^ behavior. This seems to be surprising, especially 
as far as the CMB and BAO data are concerned. But as already mentioned, for these tests 
there is just one data point and the corresponding values are negligible anyway. 
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FIG. 8: Future event horizon cutoff with ^ = 2. Left panel: one-dimensional PDF for vq. Center panel: 
one-dimensional PDF for c^. Right panel: Contour plots for the la, 2a and 3a confidence levels. 




FIG. 9: Future event horizon cutoff with ^ = 3. Left panel: one-dimensional PDF for tq. Center panel: 
one-dimensional PDF for c^. Right panel: Contour plots for the la, 2a and 3a confidence levels. 

D. Ricci-scale cutoff 

1. CPL parametrization 

According to Eq. ( 196|1 . there are four free parameters, wq, wi, tq and Hq. The minimum 
value of is used to fix the values for and Hq, such that again two free parameters are 
left. The estimations for the parameters Wq and Wi are Wq = — 1.291o'q9 "^i ~ l-15]'^oii. 
Typically, we have also 0.3 < tq < 0.5 according to the set of observational data used. The 
one- and two-dimensional PDFs are displayed in figure [TOl 
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FIG. 10: Ricci-scale cutoff with CPL parametrization. Left panel: one-dimensional PDF for wq. Center 
panel: one-dimensional PDF for wi. Right panel: Contour plots for the Icr, 2a and 3cr confidence levels. 
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FIG. 11: Ricci-scale cutoff with interaction Q = 3H(3pH- Left panel: one-dimensional PDF for rg. Center 
panel: one-dimensional PDF for (3. Right panel: Contour plots for the Icr, 2a and 3a confidence levels. 

2. Interaction Q = 3Hf3pH 

In this particular case the four free parameters are: ro, w, P and Ho (cf. Eg. |97|) . Imposing 
w = —1, and fixing Hq as described above, we obtain the following estimation for the other 
two parameters at 2a: tq = 0.351q;q3 and (3 = O.OStaoi- The corresponding PDFs are shown 
in figure [TTJ 

The different models may be properly compared among themselves and also with the 
ACDM reference model, using appropriate statistical criteria. Two options are the Akaike 
Information Criterion (AIC), and the Bayesian Information Criterion (BIC), which allow 
us to compare models with different degrees of freedom. The AIC criterion uses the formula 
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AIC = Xmin + [SQI, where z/ is the number of degrees of freedom; the BIC criterion [51 1 



is based on the expression BIC = Xmin + 2z/ In A^, where is the number of observational 
points. The smaller the resulting numbers in both expressions, the higher the quality of the 
corresponding model. The results are given in table UTTl It is obvious that the ACDM model 
is the clear winner according to both of the criteria. The worst values are obtained for the 
Ricci-scale model with CPL parametrization and for the future-event-horizon model with 
^ = 1. For the other models we find values that do not differ much from each other. 

It is convenient to classify a model with respect to the differences AAIC and ABIC 
between its AIC and BIC values, respectively, and the corresponding values for a reference 
model. This establishes a scale (Jeffreys' scale) which allows for a ranking of different models 
according to the magnitude of their differences AAIC and ABIC [52]. The smaller the 
difference to the lowest AIC or BIC values, here those of the ACDM model, the better the 
model. For differences less than 2, there is strong support for the model under consideration. 
UAAIC{ABIC) < 5 the model is still weakly supported. Models with AAIC (ABIC) > 10 
should be considered as strongly disfavored. By inspection, it follows from table IIIIl that, 
using the AIC criterion, the Hubble-scale models, and the interacting Ricci-scale model 
remain competitive with the ACDM model. The other candidates are still weakly supported. 
Applying, however, the BIC criterion, all the holographic models appear to be ruled out. 
This kind of contradiction in using different evaluation criteria is well known in the literature, 
see, e.g., |53 |. 

It should be mentioned, that the assessment of competing models also strongly depends 
on the choice of priors. Such a choice may radically change the final results from the AIC 
or BIC criteria. For example, the seven years WMAP data indicate that Hq = 71.4 for 



the ACDM model 



48j. If this result is imposed as a prior for the holographic models, all 



estimations for the AIC and BIC values change and the overall picture is modified. As an 
example, let us consider the interacting Ricci-scale model (cf. Eqs. ( ISTl) and ( !97l) ). It has 
three free parameters, Hq, tq and /3. If we fix tq = 0.4, two free parameters are left, and 
now this model has the same number of degrees of freedom as the ACDM model. Under 
this condition we obtain Xtotai — 550.99, essentially the same Xtotai ACDM model. 

Since the number of free parameters is the same, both models may be considered equally 
competitive. Hence, the imposing of priors may change substantially the final result of the 
evaluation based on the AIC and BIC selection criteria. Moreover, also priors of structure 
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Model 


FP(MP) 


SN{AIC) 


SN{BIC) 


H{AIC) 


H{BIC) 


Total(^IC) 


Total(BIC) 


ACDM 


2(2) 


546, 70 


567.99 


12.07 


18.33 


554.78 


576.16 


Hubble(J^^o = 0.25) 


3(2) 


548.59 


580.53 


14.02 


23.41 


556.62 


588.69 


Hubble(J^^o = 0.30) 


3(2) 


548.59 


580.53 


14.02 


23.41 


556.62 


588.69 


Event horizon = 1) 


3(2) 


550.49 


582.43 


14.53 


23.92 


559.03 


591.10 


Event horizon = 2) 


3(2) 


548.79 


580.72 


14.09 


23.48 


556.88 


588.95 


Event horizon(^ = 3) 


3(2) 


548.95 


580.89 


14.19 


23.58 


557.15 


589.22 


Ricci (CPL param.) 


4(2) 


551.44 


594.02 


16.01 


28.53 


559.45 


602.22 


Ricci(Q = SHPph) 


3(2) 


548.66 


580.59 


14.02 


23.41 


556.67 


588.75 



TABLE III: AIC and BIC criteria for the assessment of the considered models. In the second 

column the number of free parameter of each model is specified as well as the number of marginal- 
ized parameters (in parenthesis). The difference between the numbers of free and marginalized 
parameters leads to the number of fixed parameters. 

formation may alter the final evaluation. 

Another aspect that should be stressed is that xLin ~ seems to be a saturation value 
for any possible fitting using smooth functions. The reason is the dispersion in the SNIa and 
H{z) data. The fact that this value is grosso modo achieved by the ACDM model with just 
two free parameter seems to imply that any other model with more degrees of freedom would 
be automatically disfavored. But, using priors and other observational tests, like those from 
a perturbative analysis, may change the final conclusions again. 

VII. CONCLUSIONS 

Holographic DE models are based on a field theoretical relation between cutoffs in the 
ultraviolet and in the infrared energy regions. Different choices of the infrared cutoff give 
rise to different models. Wc have presented here a detailed analysis of models with cutoffs 
at the Hubble scale, the future event horizon and a scale proportional to the Ricci length. 
Special emphasis was put on the role of interactions with a DM component and on the 
relation between the EoS parameter and the ratio of the energy densities of DM and DE. 
We considered interaction-free limits of the dynamics and scaling solutions, i.e., solutions 
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for which the energy density ratio of the dark components remains constant. In order to 
constrain the parameter spaces for the general dynamics, a Bayesian statistical analysis has 
been performed, using data from supernovae type la, the history of the Hubble parameter, 
the position of the first acoustic peak in the CMB anisotropy spectrum and baryon acoustic 
oscillations. In detail, our results can be summarized as follows. 

i) Hubble scale cutoff: for this choice a negative EoS parameter is a pure interaction effect 
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26| . The non- interacting limit cannot represent a DE component. Independently of 
the specific form of the interaction, the energy density ratio remains constant. There is no 
direct relation between this ratio and the effective EoS parameter. For suitable interactions, 
a transition from an early matter dominated epoch to an epoch of accelerated expansion 
can be obtained . We have assumed here an interaction rate proportional to a power of the 
Hubble rate. This gives rise to the dynamics of a generalized Chaplygin gas with the ACDM 
dynamics as a special case. The results of the statistical analysis are given in subsection 
IVI Bl and in figures [5] and [6l 

(ii) Future-event-horizon cutoff: in this case, there exist a direct relation between the energy 
density ratio and the EoS parameter. On the other hand, the effective EoS parameter, i.e., 
the parameter that includes the interaction, does not depend directly on the "bare" EoS 
parameter. Assuming a power-law behavior of the energy density ratio, the cosmological 
dynamics was solved analytically for several values of the power. Using the preferred values 
of the statistical analysis in figures U\ and [9] of subsection IVI CI as well as well as in table 
Hn we found three possibilities. In the first example, the present effective EoS parameter has 
a value larger than —1. Also for a ^ 1 it will never cross the phantom divide. The second 
example has a present value of the effective EoS larger than —1 as well, but this parameter 
will become smaller than —1 for large values of the scale factor. For the third example, both 
the present and the future values of the EoS parameter are of the phantom type. Different 
from the Hubble-scale cutoff, there exists a non-interacting limit of the dynamics. Scaling 
solutions do exist as well. For all the three cutoff cases discussed here, these solutions can 
be mapped on an equivalent scalar-field dynamics with an exponential potential. All models 
with a future-event-horizon cutoff suffer from the drawback that, although they may fit the 
data, they cannot reproduce an early matter dominated epoch. 

(iii) Ricci-scale cutoff: here, there exists a linear relation between the EoS parameter and 
the energy-density ratio. The general dynamics was solved for the CPL parametrization. 
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The results of the statistical analysis of subsubsection IVID II are displayed in figure [TOl 
Furthermore, we solved the dynamics for an interaction proportional to the DE density. The 
Bayesian analysis provided us with figure [TT] of subsubsection IVI D 21 The noninteracting 
limit, which exists only for time- varying equations of state, has the interesting property that 
the energy- density ratio approaches a constant, finite value for a ^ 1 which is only about 
ten time larger than the present value. Nevertheless, the equation of state approaches that 
for dust, thus recovering an early matter dominated period. 

Tables [T] and IIIII present overviews of all the cases considered here. The results of table 
mil show, that the ACDM model remains the most favored option. According to the AIC 
criterion, the Hubble-scale models and the interacting Ricci-scale model can be considered 
as competitive since they have A AIC < 2. The other models are still weakly supported. 
However, all of the models have ABIC > 10, which indicates that they are strongly dis- 
favored according to the BIC criterion. The situation changes if certain priors are used, 
which reduce the number of free parameters to the number of free parameters of the ACDM 
model. We demonstrated, that under these conditions the relevant values for the Ricci 
case, e.g., may be competitive with that of the ACDM model. Our study has to be con- 
sidered preliminary also in the sense that it is restricted to the homogeneous and isotropic 
background dynamics. Nevertheless, we believe that a systematic investigation as performed 
in this paper, may provide an idea about basic dynamical properties of the different mod- 
els and subclasses of them. What remains to be shown for a more advanced assessment 
is to consider the perturbation dynamics and to calculate the matter power spectrum and 
the anisotropy spectrum of the CMB. For a model to be competitive one has to require 
that the preferred parameter values of the perturbative analysis coincide with those of the 
background dynamics. This issue will be the subject of a future investigation. 
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